Renormalization group approach to chiral symmetry breaking in graphene 
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We investigate the development of a gapped phase in the field theory of Dirac fermions in graphene 
with long-range Coulomb interaction. In the large- Af approximation, we show that the chiral sym- 
metry is only broken below a critical number of two-component Dirac fermions N^ = 32/7r^ , that is 
exactly half the value found in quantum electrodynamics in 2-1-1 dimensions. Adopting otherwise a 
ladder approximation, we give evidence of the existence of a critical coupling at which the anoma- 
lous dimension of the order parameter of the transition diverges. This result is consistent with the 
observation that chiral symmetry breaking may be driven by the long-range Coulomb interaction in 
the Dirac field theory, despite the divergent scaling of the Fermi velocity in the low-energy limit. 



I. INTRODUCTION 

The fabrication of single layers of carbon with atomic 
thickness has provided us with a laboratory to ex- 
plore new physics [l|, Q, as the electrons in this so- 
called graphene behave at low energies as massless 
Dirac fermions, displaying conical valence and conduc- 
tion bands 3]. Apart from its quite interesting properties 
from the applied point of view, the new material offers 
the possibility of studying an electron system that is a 
variant of quantum electrodynamics (QED) in the strong 
coupling regime, with unusual features as shown for in- 
stance in Refs. JJ-Q. 

A remarkable feature of this field theory of electrons in 
graphene is its scale- invariant character [QJ. This means 
for practical purposes that, while many-body corrections 
give rise in general to dependences on the high-energy 
cutoff, these are susceptible of being reabsorbed into the 
definition of physical quantities. Consequently, some of 
the parameters of the theory may have a nontrivial scal- 
ing in the low-energy limit. The quasiparticle weight 
is for instance renormalized, and it would be driven to 
zero if its flow were not arrested by the divergence of the 
Fermi velocity in the infrared[10j. This marginal behav- 
ior leaves anyhow an imprint in the quasiparticle decay 
rate llj , with an unconventional dependence on energy 
which has been observed experimentally [12| . 

An important phenomenon that may take place in a 
system of massless Dirac fermions is the opening of a gap 
in the regime of strong interaction. In this respect, the 
case of QED in 2-1-1 dimensions can serve as a good ex- 
ample, in which the original U{N) chiral symmetry of the 
theory with N massless two-component Dirac fermions is 
spontaneously broken below a critical number of flavors 
Nc [ij|. This chiral symmetry breaking (CSB) has been 
also studied in graphene by a number of analytical [ij- 
Il9| as well as numerical methods 20'-53] . The conclusion 
to be drawn from different approaches is that a gap can 
open up in the Dirac spectrum, though the effect may 
only appear below some critical value of N and above 
some critical interaction strength. In this picture, there 
remain however important questions to be addressed, re- 
lated to the effect of the above mentioned scaling of the 



parameters in the model. We point out in particular that 
the strength of any four-fermion interaction in the Dirac 
field theory has to be measured relative to the weight 
of the kinetic energy, that scales with the Fermi veloc- 
ity. Then, it is crucial to clarify whether the divergence 
of this parameter in the infrared may prevent the CSB 
even for a small number of Dirac fermions. 

In this paper we apply renormalization group methods 
to study the CSB in the field theory of Dirac fermions 
in graphene. We consider that this electron system is 
governed at low energies by the hamiltonian 



H = ivp / cPr Tl'i{r)j • VV^i(r) 

SttJ J ri-r2 



P(r2) 



(1) 



where {ipi} is a collection of A^/2 four-component Dirac 

spinors, tp^ = -01 70, and p(r) = Vi(r)7oV'i(r)- The ma- 
trices 7cr satisfy {7^, 7^} = 2 diag(l, —1, —1) and can be 
conveniently represented in terms of Pauli matrices as 
7o,i,2 = io'3,a3ai,a3<72) ^ (T3, where the first factor acts 
on the two sublattice components of the graphene lattice. 
Our aim is to elucidate whether a term of the type 



p„(r) =-0(r)-0(r) 



(2) 



is generated spontaneously in the hamiltonian of the elec- 
tron system. A convenient way to address this question 
is to look at the susceptibility built from that operator, 
that is, at the correlator 



n(q,w) 



-\-oo 



die^"*(rp„(q,t)p,„(-q,0)) (3) 



n(q, uj) is actually a response function measuring the re- 
action of the system under a slight difference of scalar 
potential in the two sublattices of the graphene lattice. 
A divergence of 11(0, 0) at some particular value of the 
coupling constant can be interpreted as the signal that 
Pm is getting a nonvanishing expectation value, which is 
in turn the signature of the opening of a gap in the Dirac 
spectrum. 

We will take advantage of the power of the renormaliza- 
tion group to characterize the possible singular behavior 




FIG. 1: Quantum corrections to the vertex built from tptl) 
in the large- A*' approximation, where the interaction between 
electrons is taken as the RPA dressed Coulomb potential 
(thick wavy line). 



of n as a function of e^ jvp- For this purpose, we con- 
centrate on the corrections to the vertex built from V'V'j 
as shown in Figs. [T] and |3l In the process of renormal- 
ization, /?,„ may get in general an anomalous dimension 
7^2 , modifying the naive scaling of the susceptibility. 



n(q,0)^|q 



1-27^ 



(4) 



In what follows, we apply different approaches for the 
determination of 7^2 , in order to establish the existence 
of a singular behavior in the long-distance scaling of the 
susceptibility 11. 



II. LARGE-N APPROXIMATION 

We can go beyond the usual perturbative approach in 
the coupling e^ jvp by taking formally a large number N 
of fermion flavors, to perform then the sum of all the di- 
agrams that arise to leading order in a XjN expansion. If 
we think of all possible contributions to the expectation 
value (/9m(q)?A(k -I- q)?/;'l'(k)), it is clear that the lead- 
ing corrections in 1/iV are given by the iteration of the 
exchange of electron-hole bubbles in the interaction be- 
tween the V' and '0^ fields. This amounts to adopt the 
RPA for the dressed Coulomb interaction represented in 
Fig. [Ij Introducing the polarization x(q, i^g), we get for 
the corresponding vertex function 

r(q;k) =-/o + iJ2 (27r)2 2^'" '^"''^'^^^'^° 



n=0 



Go(p + q,Wp) 



2|p-k| 



2|p-k| 



(5) 



where Gq stands for the free Dirac propagator. We recall 
that, in the case of N two-component Dirac fermions, 

X(q,c.,) = -(7V/16)qVy^4q2-c.2. 

An interesting feature of the sum in (O is that, while 
a high-energy cutoff A has to be imposed to make the 
integrals finite, all the terms show the same degree of log- 
arithmic dependence on the cutoff. By passing to imagi- 
nary frequency iujp = ujp, we can compute the divergent 



contribution to the vertex as 



7o 



n=0 ^ ^ 



cPp dUJp 



(27r)2 2tt (Zu2 -H 4p2)«/2+i 

^, V^^ 1\"„"+1 V2 ^ 2/ 



70 + 70 



79 



n=0 

arccos{g) 



d{vF\p\) 



vf\P\ 



^2Ar-yr^ 



log A 



(6) 



where g = {N/32)e'^/vp. We note that the singu- 
larity at g = 1 is only apparent, as the function in 
© can be continued analytically to g > 1 by taking 

arccos((7) — i\og{g + \/ g^ — 1). In this way, we end up 
with an expression of the vertex that becomes sensible 
for arbitrarily large values of the effective coupling. 

The divergence of the vertex F at large A has to be 
removed by absorbing the dependence on the cutoff into 



the scale Z^i of the composite field V''0|2J|- However 



this is not the only field redefinition to be accomplished, 
as the finiteness of the full Dirac propagator demands 
the introduction of a cutoff-dependent scale for the Dirac 

1 /2 

field, such that V'(A) — Zj {h)%lj^cn- The electron self- 
energy can be actually found in Ref. 'lO" to dominant 
order in the 1/iV approximation, providing the result 
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g arccos g 



TT^N 



9 VT 



log A (7) 



The cutoff-independence of the vertex F must be guar- 
anteed after multiplication by Z^ and the scale Z^2 of 
We define the renormalized vertex 
By imposing the finiteness of Fron, 



the composite field. 

as 1 rrn = Z^2Z.^i 



we obtain to leading order in the 1/A^ expansion 



1- 



n^N 



2 arccos q 
2+- ^ 



9VT 



, - logA (8) 

r 9 1 



The knowledge of Z^,2 can now be used to determine 
the anomalous scaling of the susceptibility n(q, w). This 
correlator involves two composite operators ipip, and it 
can be made cutoff-independent by multiplying each of 
them by their renormalization factor. The finite suscep- 
tibility is then nrcnlq:"^) — ■^3^2n(q, w). A renormaliza- 
tion group equation can be obtained for 11, relying on the 
independence of the susceptibility on A in the renormal- 
ized theory ;24|. We obtain from the invariance of Ilren 



^-^+^i9)§-g-^l^^]^i^^^)=0 



with the anomalous dimension 

d log Z^2 



1^^ 



-A- 



dh 



(9) 



(10) 
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FIG. 2: Phase diagram obtained to leading order in the 
\/N approximation, showing the regime with massless Dirac 
fermions (m — 0) and the phase with CSB {m j^ 0). 



and f3{g) = dg/d\ogA. In the bare theory with a 
cutoff A, it follows from dimensional analysis that the 
susceptibility 11 can be written in terms of a dimen- 
sionless function $(x) as n(q, 0) — {\q\/vF)^{vF\<l\/J^)- 
Then, neglecting in a first approximation the scaling of 
the effective coupling, the solution of ([9]) implies that 
n(q, 0) - {\q\/vF){vF\q\/A)-'^'^^^ , with the behavior an- 
ticipated in Eq. (|H) . 

The anomalous dimension obtained from ([5]) is 



Tv-^ 



TT^N 



2 arccosq 
2-\ ^ 



9VT 



r 



(11) 



and it turns out to be a monotonous, increasing function 
of g. This means that, provided that it gets sufficiently 
large, there may exist a critical value gc at which n(q, 0) 
becomes singular in the limit q — > 0. The divergence of 
this susceptibility implies a long-wavelength instability, 
which can be interpreted as the development of a nonva- 
nishing expectation value of ipip. On the other hand, the 
value of (7c depends in general on the number of flavors N. 
We can draw then a boundary marking the onset of CSB 
in (N, g) space. This line of transition, characterized by 
the condition 1 — 27^,2 = 0, is shown in Fig. [2] 

The expression (fTT|) leads to the existence of a criti- 
cal number of flavors N^, above which CSB cannot take 
place. If we let g — > 00 in that equation, we approach the 
maximum value of 7^2, from which we find N^ = 32/7r^. 
It is very suggestive that this critical N is precisely half 
the value obtained in QED in 2-1-1 dimensions |l3|. Tech- 
nically, the methods used to derive Nc in each model can- 
not be easily compared, as QED is not a scale-invariant 
field theory in that number of dimensions. On intuitive 
grounds, however, one can understand the relation be- 
tween the two values of Nc, as the photon propagating 
in QED has two different degrees of freedom. This may 
explain that twice the number of flavors are needed there 
to equally screen the interaction, in comparison to our 
model with just the scalar Coulomb potential. 





FIG. 3: Self-consistent diagrammatic equation for the vertex 
r(q;k), equivalent to the sum of ladder diagrams built from 
the iteration of the bare Coulomb interaction (thin wavy line) . 



III. LADDER APPROXIMATION 

We resort now to an approach that can better cap- 
ture the behavior of the system to the right of the phase 
diagram of Fig. [H For this purpose, it is pertinent to 
adopt a self-consistent approximation in the calculation 
of the vertex F, equivalent to the sum of ladder diagrams, 
by which the most divergent logarithmic dependences are 
taken into account at each perturbative level[25|. The ap- 
proach is encoded in the self-consistent equation shown 
in Fig. [3] The perturbative solution leads to a power 
series in the effective coupling A = e'^/8nvF, where the 
term of order A" diverges in general with the high-energy 
cutoff as log" (A) . The important point is that the set of 
diagrams considered in this way allows to implement a 



consistent renormalization of the theory, where Z^2 is 
free of nonlocal divergences, making possible a precise 
computation of the anomalous dimension 7^2 . 

A solution of the equation in Fig. [3] has been given 
in Ref. [2^ regularizing the momentum integrals with an 
infrared and a high-energy cutoff. Here, in order to fa- 
cilitate the calculation of the divergences of the vertex 
r, we define instead the field theory by analytic contin- 
uation to spatial dimension d = 2 — e. After integration 
in the frequency variable, the self-consistent equation for 
the vertex takes the form 



r(0; k) = 70 + 27rAc 



d'^p 



i2n) 



^r(o;p) 



1 



1 



IPIIP- 



(12) 



where the dimensionful coupling Aq is given in terms of 
an auxiliary momentum scale p by Aq = Ap*^ (see below) . 
Thus, powers of log A are traded by poles at e = in 
the different perturbative contributions, which are eas- 
ier to compute. In principle, Eq. p^ could also af- 
ford a nonperturbative resolution, but the computation 
of the anomalous dimension would be complicated then 
as this is obtained from the residue of the 1/e pole. Un- 
fortunately, a closed equation for that quantity cannot 
be written from Eq. (1121) . which couples the equations 
for the coefficients of the different powers of e. This is 
otherwise a natural consequence of the regularization of 
the diagrams, since the interdependence of the different 
poles is a key property of a renormalizable theory, as we 
illustrate below. 

We resort then to an iterative resolution of Eq. (|12p . 
by which we can obtain a recursion between consecutive 



orders in the power series for the vertex 

r(0;k)=7o(l + £A[,T„(k)) (13) 



n=l 



It can be easily seen that the momentum dependence of 
the different orders takes in general the form 



r„(k) 



(14) 



Inserting the n-th term of the series in the right-hand- 
side of Eq. (fT2]) , we get 



|k|("+i)^ 



27r 



27ra„ 



d'^p an 1 1 

(27r)'*|p|-p|p-k| 

d'p r(i + f) 

i27rr 0Fr(i±^) 



/•I ^-1/2(1 _2.)(-l+„e)/2 

Jo [(p-k)2x + p2(l-x)]l+W2 

1 , (4^)^/2 £(l^^) «" 



27^^ ^ r(i^) |k|("+i) 



X dx 



2;(l + (n+l)£)/2(-x „2;)(l+^)/2 



(15) 



After performing the integral in the x parameter, we find 
the relation 



fln+l =Pri+l(e) an 



with 



,„(,). i(4.)./2jXf)r(^xiMM. 



2V^ r( 



l + (rt-l)e 



)r(l-^e) 



(16) 



(17) 



In this approach, the bare vertex function can be written 
in compact form as 

00 ne ^ 

r(0;k)=7o + 7oE^"^nP^-(^) (18) 

n=l ' ' j=l 

where p is a momentum scale introduced to get the di- 
mensionless coupling X — p^'^X^. 

In the ladder approximation, it can be easily seen that 
Z^ = 1. On the other hand, the renormalization factor 
Z^2 must have the general structure 



1 + E 



c„(A) 



(19) 



The position of the different poles is determined by re- 
quiring the finiteness of Fron = Z^2T in the limit e — >■ 0. 



From the expression (jlSp . we can obtain the power series 

ci(A) = -A-log(2)A^-21og2(2)A^ 

-(flog3(2)-fiC(3))A4 

-(f log^(2) + log(2)C(3))A5 

_(2|B iog5(2) + 61og2(2)C(3) + ^C(5)) A« + . • . 
C2(A) - iA2 + log(2)A3-f flog2(2)A4 

-t-(flog3(2)-HiC(3))A^ 
+(241og^(2) + |log(2)C(3))A« + ... 
C3(A) = -i A^ - i log(2) A^ - I log2(2) A^ 

-(flog3(2)-f^C(3))A^ + ... 
C4(A) = ^A4 + ilog(2)A5 + ^log2(2)A6-f... 
C5(A) - -^A5-^log(2)A6 + ... 
C6(A) = 7^0 A6 + ... 



(20) 



and so on, with the next c„(A) starting each time with 
one more power of the coupling. 

Of all the poles, only the first can contribute to the 
anomalous dimension 7^2 . This is because the theory 
at d 7^ 2 has a finite limit A — > 00, and the cutoff only 
appears from the need to define the units of dimensionful 
quantities like p. The implicit dependence A ^ A'Aq 
leads to A(aA/aA) = eA and, following Eq. ([TOl) . 



7v= 



= -A 



dX d log Z,p-2 
dA dX 



= -A 



dci 
dA 



(21) 



In principle, the right-hand-side of (|2ip can contain con- 
tributions from higher order poles in (TTOl) . but these will 
vanish provided that dcn+i/dX = Cn{dci / dX) , identically 
for every n [27|. These are key constraints in order to 
have a renormalizable theory, since they guarantee the 
finiteness of 7^2 in the limit e — >■ 0. Quite remarkably, 
we have checked that those relations are indeed satisfied 
in our case, up to the order A^ for which we have com- 
puted the exact expression of Z^2 . 

The other important check we have made along the 
way is that Z^2 does not contain nonlocal divergences 
proportional to log(|k|/p), which appear at intermediate 
stages of the calculation. In the case of the simple pole, 
we have the result to order A® 

ci(A) = -A - log(2) A^ - 21og2(2) A^ 

-(f log3(2) + iC(3)) A^ - (f log4(2) + log(2)C(3)) A^ 
-(288 iog5(2) + 6 log2(2)C(3) + ^1.^(5)) A6 - (9f4 log6(2) 
+ f log^(2)C(3) + iC'(3) + |log(2)C(5)) A^ 



262144 
315 



log7(2) + 512 log4(2)C(3) + 2 log(2)C2(3) 



+61og2(2)C(5) + J.C(7))A«-KO(A«) 



(22) 



The general term of this series does not have a simple 
expression, but one can still obtain numerically higher 
orders of the perturbative expansion to determine the 
behavior of the function ci(A). Thus, we have computed 
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FIG. 4: Plot of the absolute value of the coefBcients c\ in 
the expansion of ci(A) as a power series in the coupling A. 



the coefficients c" of the power series in A up to order 
A^^, what is enough to establish their exponential growth 
with n. The results are displayed in Fig. |4l showing that 



ci(A) 



OO 

E 



a" A" + regular terms 



(23) 



A best fit of the asymptotic behavior at large n gives the 
value a ~ 4.5. 

The important point is the evidence that the anoma- 
lous dimension 7^2 obtained from (j2ip must have a singu- 
larity at a finite value of the effective coupling A* = 1/a. 
As one approaches this value from below, the anoma- 
lous dimension gets arbitrarily large, meaning that the 
opening of a gap is the effect that has to prevail in the 
system near A*, in spite of the upward renormalization 
of the Fermi velocity at low energies |28iJ. This applies 
in particular to the theory with small number of fla- 
vors, implying that, to the right of the phase diagram 
in Fig. [21 CSB should take place above a critical cou- 
pling e'^/8nvF = 1/a (« 0.2). 



IV. CONCLUSION 

In this paper we have applied renormalization group 
methods to analyze the development of a gapped phase 
in graphene, taking advantage of the scaling properties of 
the theory of interacting Dirac fermions in the 2D system. 



In this regard, an important effect that may question the 
breakdown of the chiral symmetry is the divergence of the 
renormalized Fermi velocity at low energies. In principle, 
the downward scaling of the effective coupling e^ /vp can 
prevent to remain above the line of the transition in Fig. 
[21 even in cases where the nominal value of the coupling 
places the model inside the region with m ^ 0. Similar 
objection for the CSB can be applied to additional lo- 
cal four-fermion interactions, as their relative strength is 
always to be measured with respect to the scale of the 
kinetic energy. One may argue however that, in a statis- 
tical formulation of the problem, there has to be a critical 
temperature for the transition to the gapped phase. The 
temperature is also a relevant scale arresting the renor- 
malization of the Fermi velocity at low energies. Then, 
it is feasible that the renormalized coupling e^/vp may 
still keep a sufficiently large value to force the transition 
at the critical energy scale. 

We note that our results in the large- A^ approach estab- 
lish that, for the physical value A^ = 4, graphene would 
remain in the gapless phase even for the largest values 
of the effective coupling attained in vacuum (e'^ /Airvp ~ 
2.2). This is in agreement with the fact that no evidence 
of transition to an insulating state has been found in 
free-standing graphene. The other important conclusion 
is that CSB must exist anyhow at sufficiently small values 
of A^, given the evidence we have obtained of a critical 
coupling at which the anomalous dimension of the order 
parameter diverges. This result could explain the obser- 
vation of a transition in Monte Carlo simulations of the 
long-range Coulomb interaction in the 2D system [20ll21|. 

A natural prediction from our analysis is that the 
gapped phase should emerge at some point in the way 
from A^ = 4 to A' = 1. The spin projection can be frozen 
for instance by applying a magnetic field, and it is ac- 
tually very appealing to think that the metal-insulator 
transition observed in that case in graphene may rest on 
this effect of CSB. It remains to be seen whether quench- 
ing also the Dirac-valley degree of freedom could lead 
to an insulating state for accessible values of e'^/vp, in 
accordance with the results of this study. 
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